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Abstract 

The evolution of the structure factor is studied during the phase-ordering 
dynamics of the kinetic Ising model with conserved order parameter. A 
preasymptotic multiscaling regime is found as in the solution of the Cahn- 
Hilliard-Cook equation, revealing that the late stage of phase-ordering is al- 
ways approached through a crossover from multiscaling to standard scaling, 
independently from the nature of the microscopic dynamics. 
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Recently we have made an extensive study |T[ of the overall time evolution in the phase- 
ordering process with scalar conserved order parameter in the framework of the Cahn- 
Hilliard-Cook (CHC) equation 0. Through the careful study of the global evolution of 
the structure factor from the very beginning of the quench down to the fully developed 
late stage, we have identified a pattern whereby the very early behavior is followed by 
an intermediate nonlinear mean-field regime before the asymptotic fully nonlinear regime 
is attained. Furthermore, the time scales of these regimes have been found to be strongly 
dependent on (i) the lengthscale considered and (ii) the parameters of the quench such as the 
amplitude A of the initial fluctuations and the final temperature Tp- The concurrence of all 
these elements gives rise to a rich and interesting variety of behaviors in the preasymptotic 
phenomenology which can be accounted for in a simple manner on the basis of a model 
introduced by Bray and Humayun (BH model) 0. 

In this paper we make an analogous study of the structure factor in the framework of the 
kinetic Ising model evolving with Kawasaki dynamics. The main point is that we observe 
also in this case the existence of an intermediate regime characterized by the mean-field 
multiscaling behavior, pointing to the generic nature of the above described structure in the 
overall time evolution. 

There is abundant experimental and numerical evidence for the universality of the late 
stage scaling in the phase-ordering process 0] yielding the following form for the structure 
factor 

C{k,t) r-. L'^{t)F{kL), (1) 

with 

L{t) ~ t^/^ (2) 

In particular the dynamic exponent z and the scaling function F{x) are independent of the 
initial condition, of the final temperature and of the nature of the dynamics, be it Langevin 
for continuous spins or Kawasaki for discrete spins. In a way this is easy to understand. 



The behavior (|T|) apphes in the large time regime, when ordered domains have formed 
and the process is dominated by the motion of interfaces, obeying an effective equation 
of motion, which presumably takes the same form independently from the nature of the 
microscopic dynamics. Now, with the results presented in this paper, we are able to extend 
the independence from the nature of the microscopic dynamics also to the preasymptotic 
regime where mean-field behavior is observed. 

In order to explain this let us summarize the results of Ref. Jll. The CHC equation of 
motion for an iV-component order parameter is given by 



dt 



dV{$) 



-V 



2, 



+ r/(x,t), (3) 



where V{(t>) = {r/2)(f? + {g/AN){(j?Y, with r < and (? > 0, while r/ is a gaussian white 
noise, with expectations 

< nix, t) > =0 

(4) 

< r/„(x, t)ri/3{x', t') > = -2TF5af3^^5{x - x')5{t - t')- 

For future reference let us denote by Mq = —r/g the square of the order parameter at 
the bottom of the local potential, which we will take 0(1). Starting from this equation, 
through a combination of the 1/A^-expansion and the gaussian auxiliary field approximation 
of Mazenko |^ , Bray and Humayun have derived a nonlinear closed equation for the structure 
factor [| 

^£^2ll = -2k^[k^ + R{t)]C{k, t) - 2^R{t)D{k, t) + 2k^TF, (5) 



where 



and 



R{t)=r + gj ^^C{k,t) (6) 



^^^' ^^ = Ij^J (^^^^ ~ ^^' ^^^^^^ " ^^' ^^^^^^' ^)- ^^^ 



{2ttYJ (27r)'^ 
Integrating formally with the initial condition 



C{k,t = 0) = A (8) 



one can write the structure factor as the sum of three pieces 



where 



C{k,t) = ACo{k,t)+TFCTAk,t) + ^CUk,t) (9) 



Co{k,t) =exp|-2fc=^ f' dt'[k^ + R{t')]\ (10) 



CTAk,t) = 2e dt'^^^^Z^ (11) 



* , Co{k,t) 
Co{k,t') 

C„i{k, t) = -2e t dt'^^^^R{t')D{k, t') (12) 

h Co{k,t') 

which are coupled together through the definitions of R and D. 

In general this is a very complicated nonlinear integral equation which must be handled 
numerically 0. However, as we have shown in Ref. [|l|, progress in the understanding of 
the solution can be made by considering the scaling properties of the three terms sepa- 
rately. Then, the first two terms, which are what is left after taking the large- iV limit, obey 
multiscaling [0. Namely one has 

Co,T,(^,t) -/:?'"" ^'Vo,T,(x) (13) 

with Ci = L{t)[kmit)L{t)]^/'^-\ Lit) = t^/^, x = kC2 and £2 = k-^{t) is the inverse of the 
peak wave vector. The two lengths £1 and £2 are related by 

^ = (dlnL)^. (14) 

The spectra of exponents a^^Xp i.^) are given by 

ao(x) = dip{x) (15) 

with 



ij{x) = 1 - {1 - x'^f (16) 

and 

2 + (d-2)ip(x) X <x* 

«T.(x) = (17) 

2 X > X* 

where x* = a/2. Finally Fq{x) = 1 while Ftp{x) = Tp/x^. 

The last term in the right hand side of (|^) has been shown by Bray and Humayun to 
obey standard scaling 

Cni{k,t) ^ L\t)Fni{x). (18) 

Summarizing, the behavior of the full solution of (H), is expected to be the outcome of the 
competition of three terms of the type 

C{k, t) ^ AL"''(^) + TpL^'^p ^""^ + {l/N)L'^ (19) 

where, for simplicity, we have omitted logarithmic factors taking Ci = L. 

As it is evident from the plot of the spectra of exponents in Fig. |1|, since ao^Tp (x) < d, 
for any finite A^ the last term in the right hand side of (|T^) eventually dominates on all 
lengthscales, leading to the asymptotic standard scaling behavior of the structure factor. 
However, due to the x dependence of the exponents, the competition takes place differently 
on different lengthscales. 

For the full analysis of the interplay of the three terms we refer to Ref. [||]. In the present 
case with d = 2, the second and third term in Eq. ([T9| ) scale in the same way and the 
interplay is only between the first and the other two together. Notice that the difference 
6a{x) = d — Q;o(x) vanishes for a; = 1, namely at the peak of the structure factor, while 
it is positive for all x 7^ 1 and in particular becomes very large for x > x*, where ao(a;) is 
negative. Therefore, even if the last two terms are bound to dominate, by modulating the 
choice of A, the relative weight of the first term may be adjusted so that multiscaling can 
be observed over a sizable preasymptotic time interval, with a spread in x about the peak 



which depends on the actual value of A. With a continuous order parameter the value of A 
can be varied at will and in Ref. ||I| it was shown that when A is very small (A -C Mq) there 
exist an observable mean-field preasymptotic behavior practically only at the peak, while 
for large values of A (A ~ Mq) preasymptotic mean-field multiscaling is clearly observed 
over the range x ^ x*. The interest of this result is that it is not just a property of the 
BH model, but the same pattern of behavior is observed in the simulation of the full CHC 
equation [|l| with A^ = 1. We now show that a very similar structure in the overall time 
evolution appears also in the dynamics of the Ising model. 

The analysis of the crossover between multiscaling and standard scaling is carried out 
by computing the effective spectrum of exponents a{x,t) that one obtains when fitting the 
structure factor to the form 

C(^,t)=£"^"'*^F(x), (20) 

with X = kC2 and C2it) = k^{t). A delicate point here is the choice of Ci. Adopting 
straightforwardly the definition of the large- A^ model given after Eq. ( ]TB| ) amounts to put 
in by hand the growth law L{t) = t^^^. On the other hand, in the large- A^ model £i is also 
given by [Co{km, t)Y^'^. Therefore we take as an unbiased choice £i = [C{km, t)Y^'^, dropping 
an inessential costant factor -F(l). In practice a{x,t) is obtained as the slope of the plot of 
lnC{x/C2(t),t) versus \nCi{t) for fixed x. When the time dependence of a{x,t) disappears 
there is scaling, which is of the standard type if a{x) is independent of x and it is of the 
multiscaling type if a{x) actually displays a dependence on x. 

We have considered an Ising system on a two dimensional lattice of size 512x512 with 
periodic boundary conditions and initially prepared in an infinite temperature configuration. 
We have then let the system evolve with Kawasaki spin-exchange dynamics after a sudden 
quench to a fixed temperature Tp below the critical temperature T^. We have computed 
the structure factor by averaging over 10^ realizations of the time histories. We have di- 
vided the entire duration of the simulation (10^ Monte Carlo steps/spin) into subsequent 
and nonoverlapping time intervals, corresponding to temporal decades, and computed the 



effective exponent in the different time intervals. 

In Fig. 1^ we have plotted the effective exponent as a function of x in the five decades 
for a quench to Tp = 0.5Tc. The curves follow the pattern of the results reported in Ref. |]I| 
for the CHC equation when A is large, as indeed it is the case for the Ising system where 
A = 1. The gross feature is that there are two markedly distinct behaviors for x < x* 
and for x > x*. To the right of x* the curves display a clear time dependence eventually 
reaching the standard scaling behavior a{x) ~ 2 in the latest time interval. To the left 
of X*, instead, the curves are bunched together displaying multiscaling with a spectrum of 
exponents which follows closely ao{x). This pattern fits quite well with the results of the 
previous analysis of the competition between the first and the other two terms on the right 
hand side of Eq. (0). The multiscaling behavior for x < x* shows that the relative size 
of the prefactors is such that, for these values of x, there exists a long time interval during 
which the structure factor is dominated by AL"°*^^''. This is the preasymptotic mean-field 
scaling which preceeds the crossover towards the eventual standard scaling |Q . The duration 
of our simulation allows to detect just the onset of the crossover and it is not long enough to 
see the definitive establishment of standard scaling. For x > x*, instead, 6a{x) = d — ao{x) 
is much too big for observing any multiscaling and the time dependence of the effective 
exponent, jointly to the fiat behavior as a function of x, shows that the crossover towards 
standard scaling is fully under way. 

The same analysis has been performed for a quench to Tp = 0.9Tc, obtaining results 
completely analogous to those of Fig. 0. Moreover we have repeated the analysis for a one 
dimensional Ising model quenched from infinite temperature to Tp = 0.5J/kB. For such 
value of Tp the equilibrium correlation length is ^(Tp) ~ 27.3 and the system orders for 
times such that the typical domain length is smaller than ^{Tp). At the longest time in our 
simulation the typical domain size, measured as the first zero of the correlation function, 
was Rg = 9.33. The behavior of the effective exponent a{x, t) in Fig. ^ is very similar to the 
one in Fig. 0, displaying a neater approach to standard scaling for x > 1. 

In summary, the study of the overall time evolution in the kinetic Ising model evolving 



with Kawasaki dynamics displays a clear crossover from multiscaling to standard scaling as 
observed in the simulation of the CHC equation with large initial fluctuations. The first 
remark is that preasymptotic multiscaling is independent from the nature of the microscopic 
dynamics, as well as asymptotic standard scaling. The second remark is that multiscaling is 
the characteristic signature of mean-field behavior. Therefore the results presented in this 
paper together with those for the CHC equation in Ref. [jll, show that a genuine feature of 
spinodal decomposition is that before reaching asymptotics the nonlinearity is governed by 
the mean-field mechanism on the large length scales. The reason of this phenomenon poses 
an interesting question for the theory of phase-ordering kinetics. It would be also interesting 
to look for systems with sufficiently large initial fluctuations in order to observe multiscaling 
experiment ally. 
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FIGURES 




FIG. 1. Spectrum of multiscaling exponents in the large- iV model with d=3, Tp = {aQ{x)) 
and Tp > (ar(x)). 
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FIG. 2. Effective exponent in the quench with Tp = 0.5Tc and d = 2. 
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FIG. 3. Effective exponent in the quencli witli Tp = O-bJ/ks and d = 1. In the time interval 
from 10 to 1000 £i = C{km,t) is almost constant. This implies that the interval of abscissae 
over which a{x,t) is determined, as the slope of lnC{x/C2{t),t) versus ln£i(t), is very narrow. 
As a consequence, the spectrum of effective exponents is very noisy and no real significance can 
attributed to the dotted curve in the plot. 
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